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Abstract 

We construct a new X-ray (2-10 keV) luminosity function of Compton-thin active galactic nuclei (AGNs) 
in the local universe, using the first MAXI/GSC source catalog surveyed in the 4-10 keV band. The sample 
consists of 37 non-blazar AGNs at z = 0.002 — 0.2, whose identification is highly (> 97%) complete. We 
confirm the trend that the fraction of absorbed AGNs with Ah > 10 22 cm~ 2 rapidly decreases against 
luminosity (L x ), from 0.73±0.25 at L x = 10 42 " 43 5 erg s" 1 to 0.12±0.09 at L x = lO 43 5 " 45 5 erg s" 1 . The 
obtained luminosity function is well fitted with a smoothly connected double power-law model whose indices 
are 71 = 0.84 (fixed) and 72 = 2.0 ±0.2 below and above the break luminosity, L* = io 43 3±0 - 4 ergs s _1 , 
respectively. While the result of the MAXI/GSC agrees well with that of HEAO-1 at L x > 10 43 ' 5 erg s" 1 , it 
gives a larger number density at the lower luminosity range. Comparison between our luminosity function 
in the 2-10 keV band and that in the 14-195 keV band obtained from the Swift/BAT survey indicates 
that the averaged broad band spectra in the 2-200 keV band should depend on luminosity, approximated 
by r ~ 1.7 for Lx ^ 10 44 erg s" 1 while T ~ 2.0 for Lx ^ 10 44 erg s _1 . This trend is confirmed by the 
correlation between the luminosities in the 2-10 keV and 14-195 keV bands in our sample. We argue that 
there is no contradiction in the luminosity functions between above and below 10 keV once this effect is 
taken into account. 
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1. INTRODUCTION 

The tight correlation between the mass of a supermas- 
sive black hole (SMBH) in a galactic center and that of the 
budge found in the local universe (Magorrian ct al. 1998; 
Ferrarese & Merritt 2000; Gebhardt et al. 2000; Marconi 
& Hunt 2003; Haring & Rix 2004; Hopkins ct al. 2007; 



Kormendy & Bender 2009; Gultekin et al. 2009) leads 
to an idea of the "co-evolution" of SMBHs and galaxies. 
Thus, understanding the growth of SMBHs is a funda- 
mental issue to elucidate the cosmological history of the 
universe in the relation with the evolution of galaxies. The 
key objects to study this are active galactic nuclei (AGNs), 
the phenomena where the SMBH gains its mass by accret- 
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ing gas. 

The most basic observational quantities to describe the 
cosmological evolution of AGNs is the luminosity func- 
tion (LF), the number density per comoving space as a 
function of luminosity and redshift. To derive the AGN 
LF, a statistically well-defined sample (unusually, a flux 
limited one) with complete identification obtained by un- 
biased surveys is required. Hard X-ray observations at 
energies above a few keV provide the most efficient and 
complete surveys to detect the whole AGNs including ob- 
scured ones (so-called "type 2" AGNs), the major class in 
this population (e.g., see Gilli ct al. 2007), thanks to its 
strong penetrating power against absorption by the sur- 
rounding material and to little contamination from stars 
in the host galaxy. In the past several years, combinations 
of hard X-ray surveys above 2 keV with different survey 
depths have revealed the evolution of LF of AGNs con- 
stituting the major part of the Cosmic X-ray background 
(CXB), which gives strong constraints on the scenario of 
the SMBH growth from z = to z ~ 5 (e.g., Ueda et al. 
2003; La Franca et al. 2005; Barger et al. 2005; Silverman 
et al. 2008; Ebrero et al. 2009; Yencho et al. 2009; Aird et 
al. 2010). 

To establish the X-ray LF of AGNs in the local universe 
is of great importance among these efforts, since it gives 
the reference for any evolution models. While Ueda et 
al. (2003), Silverman ct al. (2008), Ebrero ct al. (2009), 
and Yencho et al. (2009) have found that the luminosity 
dependent density evolution (LDDE; this term was orig- 
inally introduced by Miyaji et al. 2000) best describes of 
the X-ray AGN LF above 2 keV, Aird et al. (2010) recently 
suggest that the luminosity and density evolution (LADE) 
where the "shape" of the LF is constant over the whole 
redshift range gives a similarly good fit to their data. Ueda 
ct al. (2003) employ the local AGN LF in the 2-10 keV 
band based on the HEAO-1 all sky survey, by using the 
sample consisting of 49 AGNs compiled by Shinozaki et 
al. (2006). Sazonov & Revnivtsev (2004) 1 has also de- 
termined the local AGN LF in the 3-20 keV band from 
the RXTE/Slew survey (Revnivtsev et al. 2004), whose 
integrated volume cmissivity corrected for incompleteness 
is found to be by a factor of ^2 smaller than the HEAO-1 
result converted into the same energy band, however. 

More recently, hard X-ray surveys above 10 keV per- 
formed by the Swift and INTEGRAL satellites also de- 
termine the local AGN LF in the 14-195 keV or 15-55 
keV band (Swift; Tueller et al. 2008; Burlon et al. 2011) 
and in the 20-40 keV or 17-60 keV band (INTEGRAL; 
Bcckmann et al. 2006b; Sazonov et al. 2007), respectively. 
The advantage of these surveys is the least biases against 
heavily obscured AGNs, although the observed fraction of 
Compton thick objects with an absorption column density 
of JVh > 10 24 cm -2 is found to be as small as <5 percent in 
the total sample (Tueller et al. 2008; Burlon et al. 2011). 

1 The AGN X-ray fluxes (and consequently luminosities) used by 
Sazonov & Revnivtsev (2004) were underestimated by a factor of 
1.4 due to an error in the count rate to flux conversion (Sazonov 
et al. 2008). In this paper, we correct for this error whenever 
we refer to the results of Sazonov & Revnivtsev (2004). 



It is found the shape of the AGN LF above 10 keV as de- 
termined by Swift/BAT looks significantly different from 
the Shinozaki et al. (2006) result if the luminosity is sim- 
ply converted into the other band by assuming a typical 
AGN spectrum (characterized by a power law with a pho- 
ton index of swl.8.). The reasons of this discrepancy have 
not been understood yet. 

Thus, it is timely to revisit the local X-ray AGN LF be- 
low 10 keV from a new survey independently, in order to 
check the consistency with the previous works and solve 
the apparent contradictions among them. The Monitor of 
All-sky X-ray Image (MAXI) mission on the International 
Space Station (Matsuoka et al. 2009), currently in or- 
bit, provides a valuable opportunity for this. Hiroi et al. 
(2011) produce the first source catalog of the MAXI/Gas 
Slit Camera (GSC; Mihara et al. 2011; Sugizaki et al. 
2011) at high galactic latitudes (\b\ > 10°), by compiling 
the data in the 4-10 keV band accumulated for the first 7 
month since the start of its nominal operation. The cata- 
log contains 51 AGNs detected with a significance above 
7ct consisting of 39 Seyfert galaxies and 12 blazars. In this 
paper, we constrain the local AGN LF in the 4-10 keV 
band by using only non-blazar AGNs (i.e., Seyferts) in 
the Hiroi ct al. (2011) catalog. We also determine the in- 
trinsic distribution of absorption column density of AGNs 
(so-called Ah function) in the local universe using the 
same sample, and compare it with the previous results. 
Section 2 briefly describes the source sample and their X- 
ray spectral properties in terms of an absorption column 
density and a photon index determined from various ob- 
servatories. Section 3 describes the analysis method and 
obtained results. We discuss the implication in Section 4. 
The cosmological parameters of (Hq, Sl m , = (70hjo 
km s _1 Mpc -1 , 0.3, 0.7) arc adopted throughout the pa- 
per. The "log" symbol represents the base-10 logarithm, 
while "In" the natural logarithm. 

2. Sample 

To investigate the local LF and Ah function of AGNs, 
we collect the 37 non-blazar AGNs from the Hiroi et al. 
(2011) catalog at z — 0.002 — 0.2 that constitute a statisti- 
cally unbiased sample detected in the 4-10 keV band from 
an area of 34,000 deg 2 . Here we exclude Cen A, located at 
z < 0.002, and ESO 509-066, which has double nuclei and 
may be contaminated by nearby sources (see Table 1 in 
Hiroi et al. 2011). The four "confused" sources are ignored 
even if they contain contribution from AGNs (like NGC 
6814), which do not affect our results. As noted in Hiroi 
et al. (2011), the list of the X-ray brightest AGNs in the 
all sky has significantly changed from that of the HEAO-1 
survey performed 30 years ago; among 39 MAXI/GSC de- 
tected AGNs, only 17 objects are listed in both sample by 
Piccinotti et al. (1982) and that used by Shinozaki et al. 
(2006). The flux limit of the MAXI sample corresponds 
to 1.5 x 10~ n erg cm' 2 s' 1 (1.2 mCrab) in the 4-10 keV 
band. Figure 1 shows the log N - log S relation (integral 
form) of these AGNs in the 4-10 keV band, obtained by 
using the area curve presented in Figure 9 of Hiroi et al. 



Revisit of Local AGN X-ray Luminosity Function with MAXI 



3 



(2011). 



no 



m 

A 



0.0020 - 



0.0010 



0.0005 



0.0002 



0.0001 




Flux (ergs cm 



Fig. 1. log N - log S Relation of non-blazar AGNs in the 
4—10 keV band used in our analysis, determined from the 
first MAXI/GSC source catalog by Hiroi et al. (2011). 



and calculated luminosity from these parameters (Lx) are 
listed in the 7th to 9th columns of Table 1, respectively, 
together with the reference for the spectral parameters 
(10th column). In the conversion from the MAXI/GSC 
count rate into Lx, we consider a reflection component 
from cold, optically thick matter (Magdziarz & Zdziarski 
1995) with a solid angle of f2 = 2tt as adopted in Ueda et 
al. (2003), although this does not affect our result of the 
LF. For the remaining four targets 2 , we perform the same 
image analysis of the MAXI/GSC data in the 2-4 keV 
band as that in Hiroi et al. (2011) to obtain the hardness 
ratio between the 2-4 keV and 4-10 keV count rates. We 
first calculate the corresponding photon index T without 
considering any absorption; if we obtain T < 1.9, then we 
derive an absorption column density at the source redshift 
assuming an intrinsic power law with T = 1.9. The results 
of r and Ah with la statistical errors estimated in this 
way are also listed in Table 1 for these 4 targets. Figure 2 
shows the redshift (z) versus luminosity (Lx) plot for our 
sample. The open and filled circles correspond to those 
with a column density of log TVh < 22 and log Ah < 24, 
respectively. The optical type- 2 AGNs are further marked 
with the diagonal crosses. 



Table 1 summarizes the AGN list, where the first to 
sixth columns represent the catalog source No., MAXI 
source name, counterparts, optical type, 4-10 keV flux, 
and redshift, respectively Although it is known that us- 
ing the spectroscopic redshifts to estimate the distance of 
very nearby objects is subject to uncertainties due to the 
galaxy proper motion, we adopt these values for consis- 
tency with the analysis of the Swift/BAT AGNs in Tueller 
et al. (2008). We confirm that even if we instead adopt 
redshifts corrected for the infall into the Virgo cluster by 
Mould et al. (2000) to calculate the luminosity of AGNs 
at z < 0.01, our results of both LF and Ah function are 
little affected. Here we only distinguish between two opti- 
cal classes, "AGN1" (Seyfert 1.0-1.5) or "AGN2" (Seyfert 
1.8-2.0), for simplicity. We can regard that this AGN 
sample is nearly complete (99.3%), because 142 out of the 
total 143 X-ray sources are identified there. The flux er- 
rors due to the statistical uncertainties are better than 
14% for all the objects, and hence they are not taken into 
account in the following analysis. 

To compare our result with the previous works easily, 
we construct the AGN LF in the "intrinsic" 2-10 keV 
luminosity corrected for the absorption (i.e., before ab- 
sorption) at the source frame (hereafter represented as 
Lx)- Since we have the count rate in the 4-10 keV band 
from the MAXI/GSC survey, it is necessary to convert 
it to Lx by using the spectral information as well as 
the redshift for each source. Fortunately, we are able to 
find results of spectral fits in the 0.2-10 (or 0.5-10) keV 
band in the literature for 33 (out of 37) AGNs, which 
were obtained from data of either ASCA, XMM-Newton, 
BeppoSAX, Swift/XRT, or Suzaku. The spectral quality is 
sufficiently good in most cases, and hence we neglect their 
errors in the following analysis. The best-fit photon index 
r, absorption column density Ah (at the source frame), 
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Fig. 2. Redshift versus luminosity plot of our sample. The 
luminosity is the intrinsic (before absorption) one in the 2—10 
keV band estimated from the 4-10 keV MAXI/GSC count 
rate corrected for an absorption column density. The open 
and filled circles correspond to the X-ray unabsorbed (log TVjj 
< 22) and absorbed AGNs, respectively. Optical type 2 AGNs 
are marked with the diagonal crosses. 



3. Analysis and Results 

3.1. Analysis Method 

Our goal is to determine both the Ah function and 
absorption-corrected 2-10 keV LF of X-ray AGNs in the 
local universe. The calculation follows the same proce- 
dure as presented in Ueda et al. (2003), to which we refer 



2 IRAS 05078+1626, 2MASX J09235371-3141305, 4C +18.51, 
1RXS J213623. 1-622400 
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the reader for details. The same notation convention is 
adopted in this paper. The Ah function, f(L x , z; Nn) 
dlog Nn, represents a probability of finding an AGN with 
an absorption column density between log Ah and log Ah 
+ dlog Nn at a given luminosity, Lx, and redshift, z. For 
convenience, we assign log Nn = 20 for AGNs without 
any significant absorption, and consider only the range of 
log Nn < 24, since no Compton thick AGNs are present 
in the current sample. It is normalized as 



f 24 

/ f(L x ,z;Nn)d\ogNn = l. 

J20 



(1) 



The LF, &(Lx,z) in units of Mpc~ 3 , is defined so that 
d$(Lx, z)/d\ogLx gives the co-moving space density of 
all (Compton-thin) AGNs in a luminosity range between 
log Lx and log Lx + dlog Lx at a redshift of z. 

From the list of Nn and Lx in our sample, the best-fit 
parameters are searched for by minimizing the likelihood 
estimator defined as 



L = -2^1n 



N(N m ,L Xi , Zi ) 



JJJ N{Nn,Lx,z)d\ogNnd\ogLxdz : 



(2) 



where the suffix i denotes each object. The term 
N(Nn,Lx,z) represents the expected number from the 
survey, 

N(NmLx,z) = f(Lx,z;Nn)x 

d l {Lx ; Z) d A (zf(l + zfAz)A(Nn,Lx,z), (3) 
alogLx dz 

where d&_(z) is the angular distance, dr/dz the differential 
look back time, and A(Nn,Lx,z) the survey area, given as 
a function of flux that is calculated from Ah, Lx, and z. 
The minimization process is carried out on the MINUIT 
software package. The ler error for a single parameter can 
be estimated from the parameter range that increases the 
L value by 1. The fit applied to the unbinned data here 
cannot estimate the normalization of the LF. Hence, we 
determine it so that the expected source number agrees 
with the observed one, and estimate its relative uncer- 
tainty only from its Poisson error (l/\/A, where N = 37). 

3.2. Nn Function 

To avoid coupling between the Ah function and LF, we 
determine them step by step in the same way as Ueda ct al. 
(2003), considering the small sample size. First, we con- 
strain the Ah function by adopting the "delta-function" 
approximation for the LF only based on the sample list. 
It reduces the formula (2) to a simpler form (see equation 
(6) in Ueda et al. 2003) where the intrinsic Ah distribu- 
tion can be evaluated directly from its observed Ah his- 
togram by taking into account the Ah dependence of the 
survey area. For the four sources whose absorptions are 
estimated from the hardness ratios of the MAXI data and 
hence have non-negligible statistical errors, we take into 
account the uncertainties in Ah by introducing the "Ah 
response matrix function" as done in Ueda et al. (2003) 
(see their Section 4.1). 



As for the shape of the Ah function, we adopt a modi- 
fied version of that used in Ueda et al. (2003). The differ- 
ence from Ueda et al. (2003) is that (1) we allow such a 
case that the Ah function at log Ah < 22 is smaller than 
that at log Ah > 22, (2) we assign 4 discrete bins with the 
same width between log Ah= 20-24 for simplicity, con- 
sidering the practical difficulties to determine Ah with 
an accuracy better than dlog Ah < 1 for objects without 
good X-ray spectral data. We define the absorption frac- 
tion ip(Lx) as that of AGNs with log Ah = 22-24 among 
those with log Ah = 20-24, which is given as a function 
of luminosity. Its possible redshift dependence is ignored, 
because our sample consists of only local AGNs. The form 
of the Ah function is expressed differently for two ip(Lx) 



ranges; 



(for ip(Lx) < 



f(Lx,z;Nn) 



1 



1 - f±§ iJ>(L x ) (20 < logiVn < 21) 



and 



(for tP(Lx) > 



f{Lx,z;Nn) 



3+3e 



f - sTarV^x) 



(21 < logiV H < 22) 
(22 < logA H < 23) 
(23<logiY H <24). 



(20 < logTVH < 21) 
(21 < lo g 7V H < 22) 
(22 < logA H < 23) 
(23 <logiY H < 24). 



(4) 



(5) 



Here e defines the ratio of the Ah function in log Ah = 
23-24 to that in log A H = 22-23. It is fixed at 1.3 (in- 
stead of 1.7 as adopted in Ueda et al. 2003), according to 
the observed Ah distribution in the Swift/BAT 9-month 
survey (23/18, Tueller et al. 2008), which agrees with the 
more recent result by Bur Ion et al. (2011) . In the former 
case (equation 4), the Ah function is flat above log Ah 
= 21, while in the latter case (equation 5), the value in 
log Ah = 21-22 is taken to be the mean of those at log 
Ah = 20-21 and log A H = 22-23. The maximum absorp- 
tion fraction is ip m ax = |+f > corresponding to the case of 
/(A H ) = at log A H = 20 £ -21. 

Figure 2 clearly shows that X-ray absorbed AGNs are 
mostly found in the lower luminosity range (log Lx < 44) . 
This confirms the trend found in many previous works 
(e.g., Ueda et al. 2003; Hasinger 2008) that the absorp- 
tion fraction reduces with an increase of the AGN lumi- 
nosity. Thus, following Ueda ct al. (2003), we model the 
absorption fraction by a linear function of log Lx within 
the maximum (see above) and minimum values, which is 
taken to be 0.1. It is represented as 

tp(L x ) = min[Vw x ,max[V>44 - /?(log^x - 44), 0.1]], (6) 

where /? and ^44 are the free parameters to be determined 
through the likelihood fit. 

Table 2 summarizes the best-fit parameters of the Ah 
function and their la errors. Figure 3(a) plots the "intrin- 
sic" Ah function (corrected for the observation bias) for 
the total sample (upper), that for low luminosities of log 
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Lx < 43.5 (middle), and that for log Lx > 43.5 (lower). 
The dependence of the absorption fraction on the lumi- 
nosity is obvious. The best-fit model of the Ah function 
calculated at the mean Lx value in each region is overplot- 
ted. Figure 3(b) shows the "observed" histogram of A H 
for these 3 luminosity ranges, on which those predicted 
from the best-fit model are superposed. 

3.3. Luminosity Function 

Using the Ah function obtained above, we finally de- 
termine the local AGN LF by maximum-likelihood fit ac- 
cording to the formula (2). We adopt the smoothly con- 
nected double power law model, one of the most standard 
descriptions for X-ray AGN LFs, given as 

d$(L x ,z = 0) 



dlogLy 



A[(L X /L*yi+(L X /L*)^} 



(7) 



To implement the effect of the cosmological evolution, we 
introduce the evolution factor represented by (1 + z) pl , 



d$(Lx,z) _ d^(L x ,0) 



dlogLx 



dlogLy 



(1 + *) 



pi 



(8) 



where we fix pi = 4.2 based on the result obtained for the 
LDDE model in Ueda et al. (2003). Note that at z < 0.2 
and log Lx > 42, their LDDE model is identical with the 
pure density evolution model as represented above. 

Due to the limited sample size, we find it difficult to con- 
strain the three free parameters of the LF, 71, 72, and L* 
simultaneously. Hence, we fix the power law slope in the 
low luminosity range at three different values, 71 = 1.0, 
71 = 0.84 (the best-fit obtained from the Swift/BAT 9- 
month survey in Tueller et al. 2008), and 71 = 0.62 (the 
best-fit from the LADE model in Aird et al. 2010). The 
results of the likelihood fit for these three cases are sum- 
marized in Table 2. Figure 4 plots the best fit local LF 
determined from the MAXI survey for the case of 71 = 0.84 
(black curve). The data points are calculated by the 
Adata/Amodei method (Miyaji et al. 2001), to which the 
lcr statistical errors are attached according to the formula 
by Gehrels (1986). We find the local AGN emissivity in 
the 2-10 keV band integrated over the log Lx = 41-47 
range is W 2 -w = (1.37± 0.23) x 10 39 h 70 erg s" 1 Mpc" 3 
for 71 = 0.84. This value is close to that obtained by 
Shinozaki et al. (2006), W 2 - W = (9.3 ± 1.3) x 10 38 h 70 erg 
s _1 Mpc -3 , 3 but significantly larger than that by Sazonov 
& Revnivtsev (2004), W 2 -w = (5.2 ± 0.5) x 10 38 h 70 erg 
s _1 Mpc -3 , converted from the 3-20 keV band LF by as- 
suming a power law photon index of 1.7. 

4. Discussion 
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We revisit the local X-ray luminosity function of non- 
blazar AGNs, together with the absorption distribution 



It is larger than that presented in Section 6.1 of Shinozaki et al. 
(2006), (5.85 ±1.17) X 10 38 h 70 erg s _1 Mpc -3 . This is because 
while Shinozaki et al. (2006) calculated the volume emissivity 
using observed (i.e., absorbed) fluxes of each AGN at luminosity 
range of log Lx > 42, we here calculate it for intrinsic (de- 
absorbed) luminosities by integrating the analytical expression 
of the LF down to log Lx = 41. 



Fig. 3. (Upper) (a) Intrinsic distribution of absorption col- 
umn density (Nu function) . From the upper to lower panels, 
for the total sample, for AGNs with log Lx = 42-43.5, and for 
AGNs with log Lx = 43.5-45.5. The la error is attached to 
the data points (open circle). The lines (red) are the best-fit 
JVh function calculated at the mean luminosity in each lu- 
minosity range. (Lower) (b) The observed histograms of iVjj 
without any correction for the observational bias. The la er- 
ror (Poisson error) is attached to the data point (filled circle). 
The lines (red) show the predicted numbers in each log A'h 
bin calculated from the best-fit Nh function. 
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Fig. 4. Local X-ray AGN luminosity function in the 2-10 
keV band. The cosmological parameters of (Ho, fi m , Q\) 
= (70/170 km s _1 Mpc~\ 0.3, 0.7) arc adopted. Data 
points (black) are the observed values from the MAXI sur- 
vey based on the "N^^^/N^a^gi" estimator with error bars 
at la level. The thick solid curve (black) is the best-fit model 
with 71 = 0.84 (fixed). The thin solid curve (red) and thin 
dashed curve (cyan) are the results from the HEAO-1 survey 
by Shinozaki ct al. (2006) and from the RXTE slew survey 
by Sazonov & Revnivtsev (2004), respectively. The latter is 
converted from the 3-20 keV band to the 2-10 keV band by 
assuming a power law photon index of 1.7 for our adopted 
Hubble constant. The thick dashed (blue) and dot-dashed 
(magenta) curves are the luminosity functions derived from 
the Swift/BAT 9 month survey (Tueller et al. 2008) converted 
from the 14-195 keV band to 2-10 keV band by assuming a 
photon index of 2.0 and 1.7, respectively. 

function, based on the first source catalog of the on-going 
M AXI/GSC extragalactic survey. In spite of the fact that 
the current MAXI/GSC source sample is smaller than 
those from the HEAO-1 (49 AGNs in Shinozaki et al. 
2006) and RXTE (76 AGNs in Sazonov & Revnivtsev 
2004) all sky surveys performed in similar energy bands, 
it has some advantages to firmly establish the statistical 
properties of AGNs below 10 keV in the following points. 
(1) The sample is highly complete (99.3% = 142/143, or 
97% = 37/38 in the worst case). (2) Since we have adopted 
a relatively conservative threshold in the source selection 
in our catalog (7a), our sample is less subject to the flux 
uncertainties in the faintest end (< 14%) and is considered 
to be free from Eddington's bias as verified in log N - log 
S relation (Figure 1). This could actually be a problem 
in the sample in Shinozaki et al. (2006), who had to cor- 
rect for such biases by simulation (sec their Appendix). 
(3) The AGN fluxes are determined from the data col- 
lected from many scans (15 times per day), and hence 
can be regarded as the long-term averaged flux, less af- 
fected by short term variability than those obtained from 
a few snap-shot observations. (4) The energy band of 4-10 
keV band is more suitable for detecting obscured AGNs 
(except for Compton thick ones), thus reducing the ob- 
servation biases for the Ah function determination. It is 
expected that the MAXI/GSC 4-10 keV sample has in- 



termediate characteristics between the surveys in softer 
X-rays than 4 keV, and hard X-ray surveys in above 10 
keV. 

In fact, we find that the "observed" fraction of absorbed 
AGNs with log Ah > 22, 32% (=12/37), is higher than the 
HEAO-1 (20%=10/49) and RXTE (22%=17/76) results, 
while it is lower than that obtained from the Swift/BAT 
survey above 15 keV, 49% (=42/86; Tueller et al. 2008). 
We obtain the intrinsic Ah distribution by correcting for 
the observational biases that AGNs with heavier absorp- 
tions are harder to be detected due to the reduction of the 
count rates. The overall shape combined from both low 
and high luminosity samples (Figure 3) is well consistent 
with the Ah distribution obtained from the Swift/BAT 
survey, which show an almost flat distribution above log 
Ah > 21 with a weak peak in the log Ah = 20-21 bin 
(Tueller et al. 2008). We also confirm the strong depen- 
dence of the absorption fraction on the X-ray luminos- 
ity. We obtain the best-fit formula to describe this re- 
lation slightly different from that in Ueda et al. (2003), 
who included much fainter AGN samples in the analy- 
sis. The slope of the absorption fraction with respect to 
log L x is steep (/3 = 0.23 ±0.02 instead of (3 = 0.10 in 
Ueda et al. 2003), and reaches a higher maximum value 
(ip = 0.82 instead of ip = 0.57) at low luminosities. Such 
sharp (even sharper) change of the absorption fraction 
against the luminosity around log Lx ~43.5 is also found 
in the Swift/BAT sample (Tueller et al. 2008; Burlon et al. 
2011). The statistical uncertainty is quite large at present, 
however. It is of great importance to investigate redshift 
evolutions of the Ah function and the relation between the 
absorption fraction and luminosity by using much larger 
samples, which shall be left for future work. 

To compare with the past results obtained from sur- 
veys in similar energy bands, we overlay the best-fit 
LFs obtained by Shinozaki et al. (2006) and Sazonov 
& Revnivtsev (2004) with the thin solid (red) and thin 
dashed (cyan) curves, respectively, in Figure 4. The 
RXTE LF is converted from the 3-20 keV band into the 
2-10 keV band by assuming a photon index of 1.7 in 
our adopted cosmological parameters (Hq — 70 km s _1 
Mpc -1 instead of Hq = 75 km s _1 Mpc -1 in Sazonov & 
Revnivtsev 2004). The systematic uncertainties due to 
the choice of photon index are small in this case (18% in 
the luminosity within a range of 1.7-2.0). Here, the nor- 
malization of the LF is corrected for the maximum factor 
of incompleteness (1/0.7), assuming that the unidentified 
targets are all AGNs whose luminosity and redshift dis- 
tribution is the same as those of the identified sample. 
As already reported by Sazonov & Revnivtsev (2004), the 
RXTE result lies significantly lower than the HEAO-1 re- 
sults by a factor of ~ 2. The origin for this discrepancy is 
unclear, but we do not pursue it further in this paper. 

As clearly seen in Figure 4, our MAXI LF is closer to 
the HEAO-1 LF than the RXTE LF. In particular, it is 
in good agreement with the HEAO-1 result at high lumi- 
nosity range above log Lx — 43.5. However, the MAXI 
LF gives a larger number density at lower luminosities by 
a factor of ~2-3. By assuming a similar slope of the LF 
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in the low luminosity range (71 = 0.84- 1.0), the MAXI 
LF favors a smaller break luminosity, log L*, 42.9-43.9, 
than the best-fit HEAO-1 value (log = 44.0 ±0.4), 
though within the statistical errors. The discrepancy can 
be partially explained if the absorption fraction at these 
low luminosity range is underestimated than the reality 
in the previous work. In fact, according to the best-fit 
Ah function, the absorption fraction at log Lx = 42.5 
is estimated to be ip = 0.72 ± 0.04 in our work, while it 
is ip = 0.65 in Ueda et al. (2003). Due to the coupling 
with the Ah function in constraining the LF parameters 
through the maximum likelihood fit (see equation (2)), the 
estimated LF for all AGNs with log A H = 20-24 would 
become smaller if we assume a lower absorption fraction 
in the Ah function, because it is hard to detect objects 
with large column densities of log Ah > 23 in the 2-10 keV 
band survey, and its space density can be only constrained 
by the extrapolation from the lower column-density range. 

Comparison of our AGN LF in the 2-10 keV band 
with the hard X-ray (> 10 keV) LF determined with the 
Swift/BAT and INTEGRAL surveys provides insights on 
the broad band properties of local AGNs. Since the frac- 
tion of Compton-thick AGNs in those hard X-ray surveys 
are negligibly small, we can directly compare them with 
our result obtained for Compton-thin AGNs. In Figure 4, 
we also plot the best-fit form of the LF by Tueller et al. 
(2008) by converting the luminosity from 14-195 keV to 
2-10 keV. In this case, the assumption of the spectrum 
strongly affects the result. We adopt two photon index, 
r = 1.7 (thick dot-dashed, magenta) and T = 2.0 (thick 
dashed, blue). Obviously, the shape of the LF is not the 
same between these bands if a single spectrum index is 
assumed for all AGNs. At low luminosity range of log Lx 
^44, the normalizations of the two LFs become consistent 
with each other for T ~ 1.7, while at higher luminosities, 
the conversion with T ~ 2.0 gives a better agreement. 

This result indicates that the averaged shape of broad 
band X-ray spectra of these AGNs depends on the lu- 
minosity, in the sense that more luminous AGNs show a 
steeper slope in the 2-200 keV range on average. To con- 
firm this picture, we make the correlation plot of lumi- 
nosity between the 2-10 keV and 14-195 keV bands using 
our MAXI sample (Figure 5). Here the hard X-ray lumi- 
nosities are taken from the Swift/BAT 22-month catalog 
except for two AGNs that are not detected there; 2MASX 
J18470283-7831494 for which we refer to the Swift/BAT 
58- month catalog, and 4C +18.51, which is not yet de- 
tected by the 58 month data and has only a flux upper 
limit of 1.1 x 10~ n erg cm" 2 s _1 in the 14-195 keV band. 
We plot two lines corresponding to a power law photon 
index of T = 1.7 (solid, magenta) and T = 2.0 (dashed, 
blue). The trend that the AGNs with lower luminosities 
have flatter slope is indeed seen. 

We have shown that since the dependence of the av- 
eraged X-ray broad band spectra on luminosity makes 
the direct comparison of LFs constructed in different en- 
ergy bands (below and above 10 keV) not straightforward, 
its apparent difference in the LF shape is not "contradic- 
tion" . This effect must be taken into account when one 




42 43 44 45 



log L (2-10 keV) 

Fig. 5. Luminosity correlation between the 2-10 keV and 
14-195 keV bands, determined from the MAXI/GSC and 
Swift/BAT surveys, respectively. No K-correction and that 
for absorption is applied to the 14-195 keV luminosity. The 
open and filled circles are X-ray unabsorbed (log Nn < 22) 
and absorbed (log Nu > 22) AGNs, respectively. Optical type 
2 AGNs are marked with the diagonal crosses. The arrow de- 
notes the upper limit for an AGN that is not detected in the 
Swift/BAT 58-month survey (Baumgartncr et al. 2011). The 
dashed (blue) and solid (magenta) lines correspond to the ra- 
tio expected from a single power law spectrum with a photon 
index of 2.0 and 1.7, respectively. 

constructs a LF in a uniform way by compiling results of 
X-ray surveys performed in different energy bands. There 
are two explanations for the reasons of the luminosity de- 
pendence of the 2-200 keV spectra. Recent studies of 
nearby AGNs have suggested that the intrinsic power law 
components of Seyfert 1 are steeper than those of Seyfert 
2 (e.g., Malizia et al. 2003; Beckmann et al. 2006a; Tueller 
et al. 2008). Because of the strong dependence of absorp- 
tion fraction on the luminosity, we mostly detect only type 
1 AGNs in the surveys at the high luminosity range, lead- 
ing to the trend we see. Second possibility is the effect of 
a reflection component, which could be more significant 
at lower luminosities, as implied by the "X-ray Baldwin" 
effects (Iwasawa & Taniguchi 1993). The presence of a 
reflection hump in the spectra increases the observed flux 
in the hard X-ray band, peaked at ~ 20 keV, and hence 
the apparent slope over the 2-200 keV band becomes flat- 
ter. To distinguish these two effects, systematic studies of 
the broad X-ray band spectra of both type 1 and type 2 
AGNs at various luminosity ranges are necessary. 

5. Conclusion 

We have constructed the local AGN X-ray luminosity 
function, utilizing our new sample consisting of 37 non- 
blazar AGNs at z = 0.002-0.2 detected in the 4-10 keV 
band from the first MAXI/GSC source catalog by Hiroi et 
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al. (2011). The sample is highly complete > 97%, and is 
less subject to uncertainties in the measured fluxes com- 
pared with the past all-sky survey missions above 2 keV. 
The conclusion of our work is summarized as follows. 

• We strongly confirm the trend that there exist more 
absorbed AGNs at lower luminosities. The frac- 
tion of absorbed AGNs with log Ah = 22-24 among 
those with log Ah < 24 corrected for the observa- 
tional biases changes from 0.73±0.25 at log Lx = 
42-43.5 to 0.12±0.09 at log L x = 43.5-45.5. The 
estimated absorption distribution (Ah function) is 
consistent with the Swift/BAT and INTEGRAL re- 
sults obtained above 10 keV. 

• The shape of the intrinsic luminosity function of 
Compton thin AGNs can be fit with a smoothly 
connected double power law. For a fixed slope of 
71 = 0.84 at a lower luminosity range, we obtain a 
break luminosity of log L* = 43.3±0.4 and a higher 
luminosity slope of 72 = 2.0 ±0.2. The L* value is 
somewhat smaller than the HEAO-1 result. The in- 
tegrated emissivity over log Lx = 41-47 is found to 
be (1.37±0.23) x 10 39 erg s" 1 Mpc" 3 , which is only 
slightly larger than the previous estimate by HEAO- 
1. The space density agrees with the HEAO-1 result 
at log Lx gtsim^3.5 but is larger at the lower lumi- 
nosity range. This may be partially explained by 
the smaller biases against absorption in our survey 
in the 4-10 keV band, which lead to a better esti- 
mate of the Ah function. 

• We compare our AGN luminosity function in the 
2-10 keV band with those derived above 10 keV, 
by converting the luminosities by assuming a sin- 
gle power law spectrum. We find that the space 
densities matches with each other for r ~ 1.7 at log 
Lx < 44, while they do for T ~ 2.0 at higher lumi- 
nosities. This suggests the luminosity dependence 
of the averaged broad X-ray band spectra over the 
~ 2 — 200 keV band. The trend is indeed confirmed 
by the luminosity correlation between the MAXI 
and Swift/BAT data in our sample. 

The work is partially supported by the Ministry of 
Education, Culture, Sports, Science and Technology 
(MEXT), Grant-in-Aid No.19047001, 20041008, 
20244015, 20540237, 21340043, 21740140, 22740120, 
23000004, 23540265, and Global-COE from MEXT "The 
Next Generation of Physics, Spun from Universality and 
Emergence" and "Nanoscicncc and Quantum Physics" . 
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Table 1. Source list of AGNs used in our study with spectral and luminosity information. 
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a : observed flux in the 4-10 keV band in units of erg cm^ 2 s _i 

: intrinsic luminosity (before absporption) in the 2—10 keV band in units of erg s 
c : in units of 10 22 cm -2 



d : reference for the X-ray spectra: (1) Winter et al. (2009), (2) Risaliti et al. (2005), (3) Shinozaki et al. (2006), (4) MAXI/GSC 
hardness ratio (this work), (5) Sambruna et al. (2011), (6) Winter & Mushotzky (2010) 



Table 2. Best fit parameters of the JVh function and 2-10 keV luminosity function of AGNs in the local universe. 



A^h function 




2-10 keV Luminosity Function 




1.3 (fixed) 0.37±0.03 


0.23±0.02 


A a 7i 72 
(1.49 ± 0.25) x 10- 5 l.OO(fixed) 
(4.2 ±0.7) x 10~ 5 0.84(fixed) 2.03^;^ 
(8.8 ±1.5) x 10~ 5 0.62(fixcd) 1.95^'^ 


Lt 

43.59+°^ 
43.30t°;i 
43.03i°;i 



a : in units of [hj Mpc 3 ] 
b : in units of [h^ erg s _1 ] 
The error is ler for a single parameter. 



